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^The  article  discusses  the  problem  of  determining  the  shape  of  a turbulent  wake  /6f 
formed  behind  a self-propelled  body  In  a medium  with  a density  variable  In  the 
direction  of  action  of  gravity.  A schematic  picture  of  the  development  of  the 
wake  behind  the  moving  object  la  as  follows:  at  the  outset,  the  diffusion  Is  the 
same  In  all  directions,  and  the  wake  expands  s)nnmetrlcally;  with  Increasing  dis- 
tance from  the  object,  the  diffusion  becomes  strongly  anisotropic.  It  decreases 
in  the  vertical  direction  under  the  action  of  gravity,  and  the  wake  assumes  a flat- 
tened shape;  In  the  volume  accupled  by  the  wake,  turbulent  mixing  within  the  wake 
gives  rise  to  a more  homogeneous  density  distribution  than  In  the  ambient  medium; 
this  volume  of  liquid  Is  no  longer  In  a state  of  equilibrium  and,  acted  upon  by 
gravity,  tends  to  return  to  this  state;  the  wake  collapses, and  this  Is  accompanied 
by  Its  further  expansion  In  the  horizontal  direction  and  excitation  of  Internal  waves.  . 


The  article  deals  with  the  problem  of  the  first  stage  of  wake  development  (up 
to  collapse),  l.e.,  the  problem  of  turbulent  diffusion  In  a stratified  medium.  The 
medium  Itself  Is  at  rest.  The  description  of  the  diffusion  process  Is  based  on  the 
equations  obtained  In  Ref.  1,  which  make  It  possible  to  describe  anisotropic  diffu- 
sion. Simplifications  are  Introduced  which  are  usually  employed  In  discussions  of 
the  problem  of  turbulent  wake  propagation;  they  are  related  to  the  experimental  fact 
that  the  free  turbulent  current  zones  are  relatively  narrow.  Molecular  diffusion 
Is  neglected.  . The  density  fluctuations  are  assumed  to  be  small,  and  therefore  they 
should  be  considered  only  In  the  terms  containing  gravitational  acceleration.^ 


The  z axis  points  upward,  parallel  to  the  action  of  gravity,  and  the  x axis 
Is  parallel  to  the  motion.  In  the  body  axis  coordinate  system,  the  flow  pattern 
Is  stationary.  The  velocity  Uq  of  the  Incoming  flow  Is  assumed  to  be  much  greater 


than  the  velocity  components  in  the  wake.  The  scales  along  the  £ , i and  Z axes 

X y z 


and  the  characteristic  values  of  fluctuation  velocities  are  different.  The  above- 
mentioned  facts  make  It  possible  to  introduce  simplifications  which  are  used  in 
deriving  the  boundary  layer  equations.  In  addition,  the  following  simplifications 
are  Introduced:  (a)  third-order  moments  corresponding  to  diffusion  processes  are 
neglected  in  all  the  equations,  this  being  the  usual  assumption  in  a wake  problem; 
otherwise,  it  Is  necessary  to  Introduce  assumptions  on  third  moments  with  average 
flow  characteristics;  (b)  It  is  assumed  that  the  scale  magnitudes  of  turbulence  L and 
energy  E of  fluctuatlonal  motion  are  constant  over  the  wake  cross  section  and  vary 
only  as  a function  of  the  longitudinal  coordinate. 


Tlio  variation  of  E across  the  wake  can  be  considered  additionally. 


Tlie  system  of  equations  Is  as  follows: 
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Tlie  quantities  { } are  averaged  over  t'^e  wake  cross  section,  A = 3.86  and 

a =0.8^  and  a are  empirical  constants, 
s 

As  in  the  case  of  wake  diffusion  in  a iiomogeneous  medium, the  solution  is 
sought  in  the  form  of  exponential  dependences  on  the  longitudinal  coordinate. 

On  the  basis  of  the  system  of  equations,  the  following  estimates  can  be  obtatned 
(terms  containing  derivatives  with  respect  to  x with  multiplier  Ug  are  not  small 
quantities,  since  T ~ x,  but  their  consideration  leads  to  small  corrections); 
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The  quantity  < Uy'u/y^  — <«„'«/ )iV*,  l,e.,  in  the  approximation  considered, 
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The  following  equations  are  obtained  for  U,  E,  and  L:  ! . , 
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To  these  equations  should  be  added  the  approximate  condition^ 
^ £/y*di/dz  — const,  z- dy dz  ^ const 

which  In  the  axlsymmetrlc  case  gives 

^Ur*dr  — const 

for  a wake  with  a zero  momentum  and 

dydz const 

for  a wake  with  a constant  momentum. 

We  Introduce  the  following  relations: 


U (0,  0,  4!)/  f/o~  (x'r,  E'  -=  li  I Eu  - (x')-^  L!  U - (^')\  / J V ^ 


where  Eq,  Lq,  ^y*»  etc.  are  the  characteristic  values  of  the  quantities  and  d is 
the  diameter  of  the  body.  Tlien 


t'  t/To  - {x'T'''  ^ -^'E'  - (®') 
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Comparing  the  terms  In  the  equations,  we  have 

P -=  2 - 28,  o - 1 - fi.  V « 


At  large  distances,  the  wake  behaves  like  a two-dimensional  one,  and  this 


follows  from  the  fact  that  the  scale  t*  tends  to  a constant  value  due  to  the  rapid 


2 

dccrcace  of  <u'  > as  x'  Increases.  This  may  be  ascertained  by  Integrating  the 


equation  of  motion 
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To  find  the  solution  of  this  equation  for  specified  Initial  data  U (0,  y' , z'), 
use  may  be  made  of  the  Fourier  transformation  method^ 


U\x\  k„  k,)  - {x\  y\  z’)  exp  (-  ik,y’  --  ik,z']  dy'dt’ 
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The  equation  becomes 


^=-£/V(V  + *i*9) 


Its  solution  will  be 

X*  X* 
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Then,  using  the  method  described  In  the  paper  cited,  we  obtain  the  solution 
In  the  form 


where  ^ T S ^ [~  (^7^)  “ f /,'*)] 

ly  = 1^4  ^ e' d*'  j , /,'  = 1^4  ^ e'«p  d*'j  ^ 


The  vertical  scale  £’  tends  to  a constant  value,  passing  through  a maximum 

(by  virtue  of  the  form  of  the  function  4>(x')).  The  maximum  corresponds  to  the 
value  of  the  coordinate 


m*  = x^Nx^/d 


and  equals 
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For  a wake  with  a constant  value  of  the  momentum,  the  Initial  velocity  distri- 
bution may  be  given  in  the  form  of  the  6 function 


lU(0,y",  z")d/dz'-  £ZoD 
and  the  velocity  distribution  1s 

For  large  values  of  the  x'  coordinate,  we  obtain 
C/(.T,0, 0)  — ! V,.  E^x~\ 


In  the  case  of  a wake  with  a zero  momentum,  one  can  obtain  an  estimate  of  the 
velocity  on  the  wake  axis  for  large  x'.  Assuming  also  that  (V  )2»  1 


we  obtain 
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To  get  an  idea  of  the  pattern  of  wake  development,  we  will  use  in  the  expres- 
sions for  iy*  and  the  dependences  obtained  for  large  x.  Then 


= [2“p]‘'‘(*')‘  - (4eo)’'-(/VTo)-’'*(.VTox')’‘ 

The  distance  to  m*  covers  the  initial  stage  of  wake  development,  in  which  the 
action  of  turbulent  diffusion  may  be  considered  the  predominant  effect.  The  decrease 
in  vertical  wake  dimensions  obtained  and  the  subsequent  pattern  of  its  development 
require  further  examination.  The  solution  obtained  makes  it  possible  to  establish 
the  similarity  criterion  characterizing  the  length  of  the  initial  stage  of  wake 
development  up  to  collapse 

iVtoi*  / d " const  or  ^ ^ ^ 

1/2 

Concerning  the  factors  Lp/d  and  Uq/Dq  , they  may  be  assumed  to  have  a 

weak  dependence  on  the  Reynolds  number  and  to  be  constant  with  an  adequate  approxi- 
mation. The  constancy  of  the  ratio  L„/d  in  jets  and  wakes  is  usually  postulated  in 

^ 1/2 

the  treatment  of  these  problems.  The  constancy  of  Uq/Eq  ' over  a certain  range  of 

Reynolds  numbers  is  indicated  by  the  experimental  data  of  Ref.  6.  Thus,  the 
similarity  criterion  is  the  proportionality  of  x*/Uq,  l.e.,  of  the  time  after 
passage  of  the  object  to  the  instant  of  wake  collapse,  to  the  period  of  natural 
oscillations  of  the  medium 


This  conclusion  corresponds  to  the  experimental  data  obtained  with  models 
(Fig.  1).  Here  the- data  are  denoted  as  follows;  (1)  0 - from  Ref.  7,  (2)  O ~ 
from  Ref.  9,  and  (3)  C - from  Ref.  8. 

Comparison  with  experimental  data  gives 
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The  following  relation  Is  obtained  for  the  maximum  value  of  the  vertical 
scale  of  the  wake: 


The  wake  flattening  coefficient  Increases  with  distance  and  may  reach  a con- 
siderable value 

” = 77  f (C)dcj 

0 

The  motion  of  models  In  a stratified  medium  was  studied  In  Ref.  7 and  8.  The 
experimental  conditions  described  below  were  different,  and  the  ratios  of  the  ver- 
tical dimension  of  the  wake  to  the  diameter  of  the  model  were  found  to  be  almost 
the  same  In  these  experiments.  It  Is  obvious  from  the  solution  obtained  that  this 
agreement  corresponds  to  the  prevailing  conditions. 

The  experimental  data  according  to  Ref.  7 and  8 were: 

Uq  cm/sec  = 45,  60;  d cm  = 2.2,  15;  |dp/dz(  g cm  ^ = 0.0052,  0.001; 

T'  sec  = 0.438,  3.16;  [T'O^/d]^^^  = 3.0,  3.6 

l.e.,  the  values  of  max  I ' should  be  fairly  similar. 

z 

Comparison  with  existing  data  suggests  that  the  equations  used  make  It  possible 
to  obtain  a solution  giving  a satisfactory  pattern  of  wake  development  In  the  Ini- 
tial stage. 
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